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ABSTRACT 

It  is  well  known  that  the  Fourier  series  of  an  analytic  and  periodic  function,  truncated 
after  2N  -f  1  terms,  converges  exponentially  with  N,  even  in  the  maximum  norm.  It  is 
also  known  that  if  the  function  is  not  periodic,  the  rate  of  convergence  deteriorates;  in 
particular  there  is  no  convergence  in  the  maximum  norm,  although  the  function  is  still 
analytic.  This  is  known  as  the  Gibbs  phenomenon.  In  this  paper  we  show  that  the 
first  2N  -f  1  Fourier  coefficients  contain  enough  information  about  the  function  so  that  an 
exponentially  convergent  approximation  (in  the  maximum  norm)  can  be  constructed.  The 
proof  is  a  constructive  one  and  makes  use  of  the  Gegenbauer  polynomials  C^(a:).  It  consists 
of  the  following  two  steps. 

In  the  first  step  we  show  that  the  first  m  coefficients  of  the  Gegenbauer  expansion  (based 
on  C,J(a;),  for  0  <  n  <  m)  of  any  L2  function  can  be  obtained,  within  exponential  accuracy, 
provided  that  both  A  and  m  are  proportional  to  (but  smaller  than)  N. 

In  the  second  step  we  construct  the  Gegenbauer  expansion  based  on  C,J  ,  0  <  n  <  m 
from  the  coefficients  found  in  the  first  step.  We  show  that  this  series  converges  exponentially 
with  N,  provided  that  the  original  function  is  analytic  (thou  gh  non-periodic). 

Thus  we  prove  that  :  The  Gibbs  phenomenon  can  be  completely  overcome. 
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1  Introduction 


In  this  paper,  we  deal  with  a  prototype  of  the  Gibbs  phenomenon  and  show  how  to  eliminate 
it.  Consider  an  analytic  but  non-periodic  function  /(x)  defined  in  [—1,1].  Notice  that  /(x) 
has  a  discontinuity  at  the  boundary  x  =  ±1  if  it  is  extended  periodically  with  period  2.  The 
Fourier  coefficients  of  /(x)  are  defined  by 

m  =  y'j(x)e-‘'-'-dx.  (1.1) 

Assume  that  the  first  2N  -f  1  Fourier  coefficients  /(fc),  |fc|  <  N,  are  known  but  the 
function  /(x)  is  not.  Our  objective  is  to  recover  the  function  /(x)  for  —1  <  x  <  1  with 
exponential  accuracy  in  the  maximum  norm.  The  traditional  Fourier  partial  sum  using  the 
first  2N  -f  1  modes 


M^)  =  me"-  (1,2) 

k=-N 

does  a  poor  job;  it  produces  a  first  order  approximation  to  /(x)  with  an  error  O(^)  away 
from  the  boundary  x  =  ±1,  and  shows  0(1)  spurious  oscillations  near  the  boundary  x  =  ±1 
known  as  the  Gibbs  ph  imenon.  Thus  there  is  no  convergence  in  the  maximum  norm. 
When  one  uses  a  filter  in  the  Fourier  space, 

=  E  <TiE'/(*:)e“"  ,  (1.3) 

k=-N 

where  are  suitably  defined  real  or  complex  numbers  which  tend  to  zero  when  |fc| 

tends  to  N,  the  situation  becomes  better.  One  can  get  exponential  accuracy  away  from  the 
boundary  x  =  ±1  if  crj^  are  chosen  as  suitable  real  numbers  [7,  9,  10,  12]  or  one  can  get 
exponential  accuracy  up  to  one  boundary  x  =  — lorx  =  lif  crj^  are  chosen  as  suitable 
complex  numbers  [3].  In  these  cases  the  approximation  is  still  in  the  space  spanned  by  the 
first  2N  -t-  1  trigonometric  polynomials  and  is  a  convolution  of  the  original  Fourier  partial 
sum  with  some  filter  kernel  which  is  an  approximate  two-sided  or  one-sided  6  function;  hence 
it  cannot  be  exponentially  accurate  in  the  maximum  norm  for  —  1  <  x  <  1.  For  the  one¬ 
sided  filters  introduced  in  [3]  one  can  use  two  different  approximations  in  —  1  <  x  <  0  and 
in  0  <  X  <  1,  right-sided  for  the  former  and  left-sided  for  the  latter,  to  obtain  exponential 
convergence  globally. 

In  this  paper  we  adopt  a  different  point  of  view.  The  idea  is  the  following:  we  realize 
that  the  problem  with  the  Fourier  approximation  is  the  non-periodicity  of  the  function  and 
the  fact  that  the  functions  are  the  solutions  of  a  regular  Sturm-Liouville  problem.  In 
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[6]  it  is  shown  that  expanding  an  analytic,  non-periodic  function  f{x)  by  the  eigenfunctions 
of  a  singular  Sturm-Liouville  problem  yields  rapid  convergence.  For  example,  a  Chebyshev 
or  Legendre  expansion  of  /(x)  converges  exponentially.  Thus  if  the  first  2N  -f  1  Fourier 
coefficients  can  provide  enough  information  to  reconstruct  the  coefficients  of  an  expansion 
based  on  a  singular  S-L  problem,  we  might  recover  the  accuracy.  Unfortunately,  one  can 
not  recover  the  coefficients  of  the  Chebyshev  or  the  Legendre  expansion  within  high  enough 
accuracy. 

In  this  paper  we  show  that  from  the  first  2N  -f-  1  Fourier  coefficients  of  an  analytic  but 
non-periodic  function,  one  can  get  the  first  m  ~  iV  coefficients  in  the  Gegenbauer  expansion 
based  on  the  Gegenbauer  polynomials  C^{x),  provided  that  the  parameter  A,  appearing  in 
the  weight  function  (1  —  x^)^“2 ,  grows  with  the  number  of  Fourier  modes-  N.  We  prove  that 
this  yields  exponential  accuracy  in  the  maximum  norm. 

Our  proof  consists  of  two  separate  and  independent  steps.  The  first  step  (Section  3)  is  to 
show  that  given  the  Fourier  partial  sum  of  the  first  2N  1  Fourier  modes,  of  an  arbitrary  L2 
function  /(x),  it  is  possible  to  recover  the  partial  sum  of  the  first  m  terms  in  the  Gegenbauer 
expansion  of  the  same  function  to  exponential  accuracy  (in  the  maximum  norm)  by  letting 
the  parameter  A  and  the  number  of  terms  m  in  the  Gegenbauer  expansion  grow  linearly  with 
N.  In  this  step  /(x)  need  not  be  smooth;  any  L2  function  will  do.  We  denote  the  error 
between  the  exact  Gegenbauer  coefficients  and  those  obtained  from  the  Fourier  coefficients 
-  the  Truncation  error.  The  results  of  this  Section  are  summarized  i  n  Theorem  3.3. 

In  the  second  step  (called  the  Regularization  error),  we  prove  the  exponential  convergence, 
in  the  maximum  norm  of  the  Gegenbauer  expansion  of  an  analytic  function  when  A  grows 
linearly  with  m.  This  is  done  in  Section  4.  The  second  step  has  its  own  interest;  it  is 
an  exponential  convergence  proof  in  the  maximum  norm  for  such  Gegenbauer  expansions  of 
analytic  functions,  where  A  increases  with  the  number  of  the  terms  used  in  the  approximation. 
The  results  of  this  Section  are  summarized  in  Theorem  4.3. 

In  Section  2  we  discuss  some  results  concerning  Gegenbauer  polynomials  that  are  relevant 
to  the  proof  and  the  computations.  These  are  combined  with  the  results  of  Sections  3  and 
4  to  yield  the  main  theorem  of  Section  5  which  demonstrates  that  one  can  construct  an 
exponentially  convergent  approximation  to  an  analytic,  non-periodic  function,  from  its  first 
2N  -f  1  Fourier  coefficients.  Finally,  in  Section  6  we  demonstrate  the  theory  with  some 
numerical  examples.  Of  special  interest  is  Example  6.1,  concerning  the  function  /(x)  =  x. 
This  function  was  used  originally  (in  1898)  to  demonstrate  the  Gibbs  phenomenon. 

We  will  use  A  or  A  for  a  generic  constant  independent  of  all  the  growing  parameters 
throughout  this  paper.  The  actual  value  of  >1  or  /I  may  be  different  in  different  locations. 
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2  Preliminaries 


In  this  section  we  will  introduce  the  Gegenbauer  polynomials  and  discuss  some  of  their 
aisymptotic  behavior.  We  rely  heavily  on  the  standardization  in  Bateman  [2],  although  for 
our  purpose  a  different  scaling  might  have  been  more  natural. 

We  start  by  defining  the  Gegenbauer  polynomials  C^{x)  in  the  following  definition. 

Definition  2.1 

The  Gegenbauer  polynomial  C^(x)  is  the  polynomial  of  order  n  that  satisfies 


J^y-x^)^-^C^{x)C^{x)dx  =  0  k^n 


and  (for  A  >  0) 


_  r(n  +  2A) 
n!r(2A)  ■ 


Note  that  the  Gegenbauer  polynomials  thus  defined  are  not  orthonormal.  In  fact,  the 
norm  of  C^(x)  is  given  by  the  following  lemma. 

Lemma  2.1 

The  Gegenbauer  polynomials  defined  above  satisfy 


£(1  -  x=)'-ici(x)C»(x)rfi  = 


where 


r(A  +  |) 
r(A)(n  +  A) 


For  the  proof  see  [2,  page  174]. 


We  are  ready  now  to  deal  with  the  asymptotics  of  the  Gegenbauer  polynomials  for  large 
n  and  A.  For  this  we  need  Lemma  2.2  and  2.3. 


Lemma  2.2  (Stirling) 

For  any  number  i  such  that  i  >  1  we  have 


r(x  +  1)  <  (27r)^x^+^e-*eA 

r(x  +  l)  >  (2jr)U^+ie-=' 
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(2.5) 

(2.6) 


□ 


Lemma  2.3 

There  exists  a  constant  A  independent  of  A  and  n  such  that 


(2.7) 

(2.8) 

The  proof  follows  from  (2.4)  and  the  Stirling’s  formula  (2.5)-(2.6). 


□ 


Finally  we  would  like  to  quote  the  Rodrigues’  formula  [2,  page  175]. 


Lemma  2.4 

The  Gegenbauer  polynomials  are  explicitly  given  by 


(1  -  n)£;  [(1  - 

where  G(A,n)  is  defined  by 


G(A,n) 


r(A -k  |)r(n  +  2A) 
r(2A)r(n  +  A  +  i)  ■ 


(2.9) 


(2.10) 

□ 


3  The  Truncation  Error 

In  this  section  we  consider  an  arbitrary  L2  function  /(x)  defined  in  [—1,1].  We  assume  that 
the  first  2N  +  1  Fourier  coefficients  f{k),  as  defined  in  (1.1),  are  given.  We  are  interested 
in  recovering  (within  exponential  accuracy  in  the  maximum  norm)  the  first  m  coefficients  in 
the  Gegenbauer  expansion  of  /(x), 

00 

m  =  'Ef'wci'{x),  (3.1) 

1=0 

where  the  Gegenbauer  coefficients  are  defined  by 
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(3.2) 


/■'(')  =  A  Ijl  -  f(x)C!-lx)dx 

with  hf'  given  by  (2.4). 

Since  we  do  not  know  the  function  /(x),  but  rather  its  truncated  Fourier  series  /Ar(x) 
(defined  in  (1.2)),  we  have  only  an  approximation  to  f^{l)  which  we  denote  by  given 

by 


s"(0  =  A Iji  -  Mx)C>-{x)dx 


(3.3) 


Notice  that  g^{l)  depends  on  N.  At  this  stage  we  would  like  to  define  the  truncation  error 


TE{X,m,N) 


=  max 

-l<r<l 


D/'(0  -  s\l))C^(x) 


1=0 


(3.4) 


The  truncation  error  is  the  difference  between  the  Gegenbauer  expansion  (with  m  terms) 
of  the  function  /(x)  and  that  of  the  truncated  Fourier  series  It  measures  the  error  in 

the  finite  Gegenbauer  expansion  due  to  truncating  the  Fourier  series. 

In  the  next  two  theorems  we  bound  the  truncation  error  in  terms  of  N ^  the  number  of 
given  Fourier  coefficients,  m  the  number  of  Gegenbauer  polynomials,  and  A. 


Theorem  3.1 

If  /(x)  is  an  L2  function  on  [—1,1],  then  there  exists  a  constant  A  which  is  independent 
of  A,  m  and  N ^  such  that  the  truncation  error  defined  in  (3.4)  satisfies  the  following  estimate; 


TE{X,m,N) 


^  (m  +  A)r(m  +  2A)r(A) 
-  (m  -  l)!r(2A) 


(3.5) 


Proof: 

As  a  first  step  we  consider  the  special  function  /(x)  =  with  |n|  >  N.  In  this  special 
case  //v(x)  =  0  and  we  obtain 

(f  (0  -  j"(0)C/(l)  =  J\l  -  x^)^-ie"”Ct(x)dx  .  (3.6) 

Roughly  speaking,  one  can  argue  that  this  integral  is  rapidly  decreasing  when  n  increases. 
This  results  from  the  fact  that  the  integral  is  proportional  to  the  n’th  Fourier  coefficient  of 
the  function  (1  -  x^)^“^C/^(x)  which  is  analytic  and  heis  A  periodic  derivatives.  It  is  nice  to 
know  that  an  explicit  expression  of  this  integral  appears  in  the  literature  [2,  page  213): 
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4 /',(>  I'(A)  i‘{l  +  (3.7) 

wli«ro  J^{x)  i«  tho  De««cl  function.  Since  \Ju{!ii)\  ^  1  for  all  x  and  ^  0  [1,  page  302],  we 
have,  for  0  ,n, 


whore  in  the  second  step  we  used  the  formula  (2.2)  for  and  in  the  last  step  we  used 

the  fact  that  an  increasing  function  of  /. 

We  now  return  to  the  general  function  f{x),  which  satisfies 


M  -  /six)  =  2  . 

Since  J{x)  is  an  function,  its  Fourier  coefficients  }{n)  arc  uniformly  bounded, 


l/(n)|  ^  A  . 


We  thus  have,  using  the  result  for  the  special  case  in  (3,8), 


< 


(m -h  A)r(m -f  2A)r(A)  ^ 

^Jr(2A)  ini'yvUH/ 

r(m  +  A)r(m  +  2A)r(A)  /  2 
m!r(2A)  Kitn) 


for  all  0  <  /  <  w. 

We  can  now  estimate  the  truncation  error  (3,4)  by 


TEiX,rn,N)  ^ 
< 

< 


.„.,(m  +  A)r(Tn  +  2A)r(A)  / 

^  mir(2A)  [wn) 


(3.10) 


(3.11) 
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where  In  tlie  Hecotid  step  we  uued  the  f<tct  that  |C'/'(aj)|  ^  C^{i)  for  all  -1  <  a;  ^  I 
[2,  page  20(1],  and  in  the  third  itep  we  it»cd  (3.11), 

The  theorem  Ib  now  proven. 


□ 


Remark  3.X  The  approximate  Oegenbaiier  coefficients  ^^{l)  in  (3,3)  can  he  explicitly 
expresBcd  in  terms  of  the  Fourier  coefficients  f{k)  as 

J‘(0  =  W(0)  +  r(A)<'((  +  A)  ■£  f{k)  .  (3.12) 

Equation  (3,12)  follows  immediately  from  the  definition  (3,3)  of  ^^{l)  and  the  integration 
formula  (3,7), 


□ 

For  fixed  A,  the  truncation  error  (3,6)  decays  algebraically  as  0  However,  if  both 

A  and  rn  grow  linearly  with  N,  the  truncation  error  can  be  made  exponentially  small.  This 
is  stated  in  the  following  theorem. 

Thgorgm..3i2 

If  A  =s  aN  and  m  «=  f3N  where  a  and  /5  arc  positive  constants,  then  the  truncation  error 
defined  in  (3.4)  satisfies 


TE{aNJN,N)  ^  ANY 

where 

^  **  (27rc)®a®/3F^ 

In  particular,  if  a  ■  then 

q  s=  «  0.8  <  1  . 


(3,13) 


(3.14) 


EiftfiCi 

We  use  Stirling’s  formula  (2,6)'(2.6)  to  obtain,  from  the  previous  estimates  on  the  trun¬ 
cation  error  in  (3.6)  and  some  simple  algebra. 
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Wc  would  like  to  point  out  that  we  choose  a  =  in  Theorem  3.2  simply  to  show  that  ii 
is  possible  to  obtain  exponentially  small  truncation  errors,  This  may  not  be  the  best  choice 
in  practice.  We  can  easily  verify  that,  for  fixed  a,  (3,14)  defines  a  q  which  is  an  increasing 
function  of  This  is  not  surprising  since  the  truncation  error  should  be  bigger  if  there  are 
more  terms  in  the  Gegenbauer  expansion  to  approximate.  However  wo  will  see  in  the  next 
section  that  the  regularization  error  will  be  smaller  if  m  is  bigger.  In  practice  one  might  try 
to  choose  m  to  maintain  some  balance  between  these  two  errors.  For  a  fixed  /?,  (3.14)  defines 
a  q  which  attains  its  minimum  at 

ir -20  + 
a  » . 

\i  0  ^  For  example,  if  0  a  and  a  »  1,33  is  chosen  according  to  (3.16),  then  q  given  by 
(3,14)  is  approximately  0,49,  much  smaller  than  the  minimum  value  0,8  obtained  with  the 
restriction  a  0. 

Wc  summarize  the  results  of  this  Section  in  the  next  theorem. 


Theorem  3.3  The  exponential  decay  of  the  truncation  error 

bet  f{x)  be  an  /vaf-l,  1]  function,  and  f{k),  -N  ^  k  <  N,  its  Fourier  coefficients  defined 
in  (1.1).  Let  /''(/)  be  the  Gegenbauer  expansion  coefficients  of  f{x)  defined  in  (3,2),  and  let 


be  the  Gegenbauer  coefficients  of  the  truncated  Fourier  series  /six)  defined  in  (1,2) 
with  arc  given  explicitly  in  (3.12). 

Then;  if  A  =  rn  =  0N,  where  0  <  the  truncation  error  decays  exponentially  with 
the  number  of  Fourier  modes  N,  i.e. 


rn 

TFA0N,mN)  =  I  E(/*(0  -  ^\0)C/‘(*)l  S  (3.16) 
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with 


□ 

4  The  Regularization  Error 

In  this  section  we  would  like  to  establish  error  estimates  for  approximating  an  analytic 
function  f{x)  on  [—1, 1]  by  its  Gegenbauer  expansion  based  on  the  Gegenbauer  polynomials 
C^(a:).  Since  our  goal  is  to  remove  the  Gibbs  phenomenon  we  will  use  the  maximum  norm. 
In  the  last  section  we  have  shown  that  we  can  get  the  Gegenbauer  partial  sum  of  the  first 
m  terms  of  any  L2  function  from  its  Fourier  partial  sum  of  the  first  2N  +  1  modes  with 
exponential  accuracy  in  the  maximum  norm,  if  A  and  m  are  both  growing  linearly  with  N . 
Thus  in  this  section  we  will  consider  the  case  of  large  A  and  m. 

We  will  assume  that  f{x)  is  an  analytic  function  on  [—1,1]  satisfying  the  following  as¬ 
sumption. 


Assumption  4.1 

There  exists  constants  p  >  1  and  C(p)  such  that,  for  every  >  0, 


<C(f>)^.  (4.1) 

This  is  a  standard  assumption  for  analytic  functions.  Here  p  is  actually  the  distance  from 
[—1,1]  to  the  nearest  singularity  of  f{x)  in  the  complex  plane  (see  for  example  [8]).  The 
assumption  can  be  modified  using  the  techniques  m  [5]. 

Let  us  consider  the  Gegenbauer  partial  sum  of  the  first  m  terms  for  the  function  f{x) 

A 

given  by  (3.1),  with  the  Gegenbauer  coefficients  /''(/)  defined  by  (3.2).  We  want  to  estimate 
the  regularization  error  in  the  maximum  norm, 


max 


dx'^ 


(i) 


RE{\,  m)  =  max 


1=0 


We  start  by  estimating  the  Gegenbauer  coefficients  /^(/). 


(4.2) 


Lemma  4.1 

The  '^egenbauer  coefficient  as  defined  in  (3.2),  of  an  analytic  function  satisfying 

the  Assumption  (4.1),  can  be  bounded  by 
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A  Cip)r{x  +  \)Vii  +  2X) 
/i^(2/))'r(2A)r(/ +  A  + 1)  ■ 


(^.3) 


Proof; 

We  start  by  using  the  definition  (3.2)  for  /'^(/).  Wc  replace  the  term  (1  —  x'^)^~^Ct'{x) 
by  the  Rodrigues’  formula  (2.9)-(2.10)  to  get 

where  C'(A,/)  is  defined  in  (2.10).  Integrating  by  parts  /  times  we  get 

Wc  now  use  the  Assumption  (4.1)  to  estimate  the  derivative  ^{x),  thus  obtaining 

1^(01  S  /' (1  -  x^^-idx  . 

Since  Co  (a:)  =  1,  the  remaining  integral  is  simply  and  can  be  obtained  from  (2.3)'(2.4): 

,  ;A.n,  .  GiXJ)C{p)y/^T{U\  +  \) 

_  4A2y(/  +  A)r(/  +  A) 

and  finally  using  the  definition  of  C(A,/)  from  (2.10)  wc  get  (4.3). 


□ 


The  estimate  (4.3)  can  be  used  naturally  to  get  an  estimate  in  the  weighted  L2  norm. 
This  would  have  been  more  transparent  if  wc  had  adopted  an  orthonormal  Gegenbauer  basis 
rather  than  (2.2).  However  the  aim  of  this  paper  is  to  establish  estimates  in  the  maximum 
norm,  This  is  stated  in  Theorem  4,1. 


Theorem  4tl 

If  f{x)  is  an  analytic  function  on  (-1, 1]  satisfying  the  Assumption  (4.1),  then  the  regu¬ 
larization  error  defined  in  (4.2)  can  be  bounded  by 


RE{X,m)  <  A 


C{p)r{X  +  \)T{m  +  2X-\-l) 
mv/A(2/))”*r(2A)r(m  +  A)  ' 


(4.4; 


E£0.0Il 

By  (4,3)  and  (2.8)  wc  obtain 
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(4.5) 


About  the  size  of  the  regularization  error  (4.4)  when  A  depends  linearly  on  in,  as  is  the 
case  in  Theorem  3.2  for  truncation  errors,  we  can  state  the  following. 

Theorem  4.2 

If  A  =  ')7f)  where  7  is  a  positive  constant,  then  the  regularization  error  defined  in  (4.2) 
satisfies 

RE(-im,rn)  <  Aq^  (4.6) 

where  q  is  given  by 
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_  (1  +  27)^+^^ 

^  ~  /921+2')'7T'(1  +7)i+t' 

which  is  always  less  than  1.  In  particular,  if  7  =  1  and  m 
constant,  then 

RE{(5N,^N)<  Aq^  (4.8) 

with 


Proof; 

VVe  use  Stirling’s  formula  (2.5)-(2.6)  to  replace  the  Gamma  functions  in  (4.4).  A  little  bit 
of  algebra  brings  us  to  (4.6)-(4.7)  if  A  =  7m.  Notice  that  the  constant  A  in  (4.6)  contains 
the  contribution  of  p  related  terms.  It  is  easy  to  verify  that  q  defined  by  (4.7)  is  an  strictly 
increasing  function  of  7  and  tends  to  A  <  1  when  7  tends  to  infinity.  Hence  we  have  9  <  ^  <  1 
for  all  7  >  0.  As  for  the  proof  of  (4.8)-(4.9),  we  simply  plug  in  7  =  1  and  m  =  into 
(4.6)-(4.7). 

□ 


Finally  we  summarize  the  results  of  this  Section  in  the  following  theorem. 

Theorem  4.3  The  exponential  decay  of  the  regularization  error 

A  . 

Let  f{x)  be  an  analytic  function  on  [—1, 1]  satisfying  assumption  (4.1).  Let  /''(/),  0  <  /  < 
m  be  its  Gegenbauer  coefficients  defined  in  (3.2).  Furthermore  assume  that  X  —  m  =  0N. 


Then 


max  f{x)-'£pO)Ct{x)  <4,"'  (4.10) 


where 


□ 
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5  The  Main  Theorem 


In  this  Section  we  bring  the  main  theorem  demonstrating  that  one  can  construct  an  ex¬ 
ponentially  convergent  (in  the  maximum  norm)  approximation  to  an  analytic,  non-periodic 
function,  from  its  first  2N  +  1  Fourier  coefficients.  The  method  is  indicated  in  the  last  two 
sections.  Namely,  first  we  get,  from  the  Fourier  coefficients,  an  approximation  to  the  first 
m  =  j3N  Gegenbauer  coefficients.  This  approximation,  by  virtue  of  the  discussion  in  Sec¬ 
tion  3,  converges  exponentially  f  ast  to  the  true  coefficients,  provided  that  A  grows  with  N . 
We  then  construct  the  partial  Gegenbauer  expansion  which  converges  exponentially  to  /(i) 
by  virtue  of  the  discussion  in  Section  4. 

We  are  ready  to  state  our  main  theorem. 

Theorem  5.1  Removal  of  the  Gibbs  Phenomenon 

Consider  an  analytic  and  non-periodic  function  f{x)  on  [—1,1],  satisfying 

Ak  f  jUI 

-KS.  -  ^<"^7  "  -  ‘ 

Assume  that  the  Fourier  coefficients 

f{k)  =  i  £ 

are  known  (or  —N  <  k  <  N. 

Let  ^^(/),  0  <  /  <  m  be  the  Gegenbauer  expansion  coefficients  of  /^{x)  = 
explicitly  given  by 

j"(0  =  W(0)  +  r(A)i'((  +  A)  ■£  m  .  (5.2) 

0<\k]<N 

Then  for  A  =  m  =  where  /?  <  we  have 


max^  /(x)  -  EsWM  <  AN^It  + 


Proof; 

We  start  by  noting  that  formula  (5.2),  for  the  approximate  Gegenbauer  coefficients  g^{l), 
follows  from  Remark  3.1. 

In  order  to  establish  (5.3),  we  introduce  the  Gegenbauer  coefficients  of  the  function  f{x) 
and  denote  them  by  they  are  defined  in  (3.2). 
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We  have 


max 

-1<X<1 


m  j 

1=0  I 


<  max 

“  -1<X<1 


max 


/=0 


+ 


1=0  1=0 


The  first  term  is  the  regularization  error  and  has  been  estimated  in  Theorem  4.3.  The 
second  term  is  the  truncation  error  and  has  been  estimated  in  Theorem  3.3.  The  theorem  is 
thus  proved. 


□ 


A  few  remarks  are  in  order: 

Remark  5.1  The  proof  is  constructive.  Given  2N  +  1  fourier  coefficients  f{k), 
—  N<k<N,  one  gets  explicitly  the  approximate  Gegenbauer  coefficients  g^{l)  ,  0  <  I  <  m 
and  the  Gegenbauer  series  can  be  explicitly  constructed. 

Remark  5.2  The  reconstruction  is  not  optimal,  since  no  effort  has  been  spent  to  optimize 
the  parameters. 

6  Numerical  Results 

In  this  section  we  demonstrate  the  theory  using  numerical  examples.  We  implement  the 
method  in  the  following  way.  Assume  that  the  first  2N  +  1  Fourier  coefficients  of  /(x),  f{k) 
for  —  jV  <  k  <  N,  as  defined  in  (1.1),  are  given.  We  compute  the  approximate  Gegenbauer 
coefficients  for  0  <  /  <  m,  defined  in  (3.3),  using  the  following  formula  given  in  Remark 
3.1: 


j*(0  =  W(0)  +  r{A)!'(/ +  A)  f{k).  (6.1) 

0<|fc|<iV  VTTfc/ 

We  compute  the  Bessel  function  Ju{x)  using  an  IMSL  routine.  Once  the  approximate 
Gegenbauer  coefficients  g^{l)  are  obtained,  we  can  compute  the  approximation  to  f{x)  by 
directly  summing 


siW  =  E9'(')C,V)  (6.2) 

1=0 

as  long  as  we  can  compute  the  Gegenbauer  polynomial  Ci'{x)  accurately.  The  formula  we 
used  to  compute  C^{x)  is 
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(6.3) 


I 


C,^cosO)=S 

k=0 


r{k  +  A)  ru  -k  +  x) 

-  —  cos 

A:!r(A)  (/-jt)!r(A) 


n-2k)0 


which  can  be  found  in  [1,  page  175]. 

We  remark  that  the  implementation  described  above  are  subject  to  round  off  effects  for 
large  A  and  m.  We  use  a  Cray-YMP  to  carry  out  all  the  computations.  Our  implementation 
can  give  accurate  result  only  when  the  error  is  no  smaller  than  10“®.  And  we  will  show  results 
only  in  those  cases.  A  better  way  to  implement  this  method  might  be  through  Chebyshev 
polynomials. 

For  the  purpose  of  testing  we  consider  the  following  examples. 


Example  6.1  f{x)  =  x 

This  is  the  original  example  used  in  1898  to  show  the  famous  Gibbs  phenomenon.  For 
this  particular  function  there  is  no  regularization  error  as  long  as  m  >  1,  hence  all  the  errors 
observed  result  from  the  truncation  error. 

In  Figure  1,  we  show  the  exact  function  f{x)  =  a:  in  the  solid  line,  the  Fourier  partial  sum 
fj\/{x)  with  N  =  4  in  the  short  dashed  line,  and  the  approximation  through  the  Gegenbauer 
polynomials  as  defined  in  (6.2),  with  N  =  A  and  m  =  A  =  jn  the  long  dashed 

line.  We  can  clearly  see  that  the  Fourier  partial  sum  fN{x)  shows  the  Gibbs  oscillations, 
while  the  approximation  g^{x)  through  the  Gegenbauer  polynomials  is  uniformly  accurate. 
In  Figure  2  we  draw  the  errors  in  the  maximum  norm,  with  a  logarithmic  scale,  versus  N 
(recall  that  2N  +  1  is  the  number  of  Fourier  modes  given),  both  with  m  —  \  =  \N  (squares) 
and  with  m  =  and  A  determined  by  (3.15)  (circles).  This  picture  confirms  our  estimates 
in  (3.13)  through  (3.15)  for  the  exponential  convergence  of  the  truncation  errors. 

Next,  we  show  the  point-wise  errors  of  the  approximations  with  m  =  X  =  ^iV  for 
N  =  4,8,16,32  in  Figure  3,  and  the  point-wise  errors  with  m  =  X/g  and  A  determined 
by  (3.15),  for  =  4,8, 16,  in  Figure  4,  again  both  in  logarithmic  scales.  We  can  observe  ex¬ 
ponential  convergence  both  inside  the  interval  and  near  the  boundary,  although  the  absolute 
error  is  several  magnitudes  smaller  inside  the  interval  than  at  the  boundary. 

We  have  also  run  the  test  for  /(x)  =  x"  with  n  =  2  through  9,  obtaining  similar  results 
(not  shown  here). 

Example  6.2  /(x)  =  cos[1.47r(x  —  1)]. 

This  function  satisfies  the  Assumption  (4.1)  with  arbitrary  p  >  1.  Both  truncation 
error  and  regularization  error  will  be  present.  The  choice  m  =  X  =  Xjg  gives  exponential 
convergence  as  expected.  The  choice  m  =  A  =  gives  better  results  in  this  case. 
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In  Figure  5  we  draw  the  errors  in  the  maximum  norm  versus  N,  with  a  logarithmic  scale, 
both  with  m  =  X  =  (squares)  and  with  in  =  A  =  (circles).  The  choice  of  m  = 
and  A  determined  by  (3.15)  does  not  converge  in  this  case,  indicating  that  (3.15)  is  probably 
a  bad  choice  for  the  Regularization  errors. 

Finally,  again  in  logarithmic  scales,  we  show  the  point-wise  errors  of  the  approximations 
with  m  —  X  =  for  N  =  24,36,52,76  in  Figure  6,  and  the  point-wise  errors  with  m  = 
A  =  for  N  =  16,24,36,52  in  Figure  7.  We  again  observe  exponential  convergence  both 
inside  the  interval  and  at  the  boundary,  and  the  several  magnitudes  difference  in  the  absolute 
errors  inside  the  interval  and  at  the  boundary. 


7  Concluding  Remarks 

We  considered  in  this  paper  a  prototype  of  the  Gibbs  phenomenon  and  demonstrated  that 
this  effect  may  be  overcome.  We  show  that  if  f{x)  is  an  analytic  but  non-periodic  function 
then  the  knowledge  of  its  first  27V  +  1  Fourier  coefficients  provides  accurate  enough  infor¬ 
mation  so  that  the  first  m  ~  TV  coefficients  in  the  expansion  of  f{x)  by  the  Gegenbauer 
polynomials  Cj^{x)  (A  ~  TV)  may  be  obtained.  In  equation  (5.2)  the  formula  that  relates  the 
Gegenbauer  coefficients  to  the  Fourier  coefficients  is  presented. 

There  are  several  topics  which  are  not  addressed  in  this  paper  and  will  be  discussed  in  a 
future  paper: 

(1)  Other  choices  of  polynomials  such  as  the  Laguerre  polynomials  and  Hermite  polyno¬ 
mials; 

(2)  Algebraic  convergence  for  but  not  analytic  functions; 

(3)  The  case  with  exponential  recovery  in  a  sub-interval  [a, 6]  of  [—1, 1].  That  is,  given 
the  first  2N  -f  1  Fourier  modes  for  a  function  defined  on  [-1,1],  find  its  Gegenbauer  partial 
sum  of  the  first  m  terms  based  on  the  scaled  Gegenbauer  polynomials  on  [a,b],  in  which 
the  function  f{x)  is  analytic.  This  would  allow  one  to  handle  multiple  discontinuities  and 
discontinuities  of  unknown  location; 

(4)  The  procedure  realized  in  the  Fourier  space  as  a  convolution,  similar  to  the  approach 
used  in  [3]; 

(5)  Efficient  and  stable  numerical  implementation  of  those  methods; 

(6)  Extension  to  Collocation  techniques  and  other  than  the  Fourier  basis  functions. 
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Figure  1:  The  function  f{x)  =  x  (background  Holid  line);  tlie  Fourier  partial  sum  /six) 
defined  in  (1.2)  with  N  <=  A  (short  dashed  line),  and  the  approximation  defined  in 
(6.2)through  the  Gegenbauer  polynomials  with  N  =  A  and  m  =  X  =  \N  in  the  long  dashed 


Figure  2:  For  the  function  f{x)  m  x,  the  errom  in  the  maximum  norm,  max-i^v^i  \fix)  - 
gm{x)\  wiiero  gm{x)  is  defined  in  (6.2),  in  a  logarithmic  scale,  versus  N  (2A^  +  1  is  tlie  number 
of  Fourier  modes  given).  The  squares  arc  for  the  case  with  m  ^  X  sz  and  the  circles  are 
for  the  case  with  mss  \N  and  X  determined  by  (3.15). 
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Logarithmic  Errors  in  Maximum  Norm 


Figure  6!  For  the  function  /(»)  ■  co»[l  An{cc  -I)],  the  error*  in  the  maximum  norm, 
max.i<£«iji !/(»)  -  whore  is  defined  in  (6,2),  in  a  logarithmic  *calc,  vcr»ui 

N  {2N  1  i*  the  number  of  Fourier  mode*  given),  The  square*  are  for  the  ca#c  with 

m  ■>  A  ■  and  the  circles  are  for  the  case  with  m  m  X  m 


Logarithmic  Pointwise  Errors 


Figure  7:  For  the  function  f{x)  =*  co8[1.4ff(a;— 1)],  in  a  logarithmic  scale,  the  point-wise  errors 
l/(®)  where  gm{x)  is  defined  in  (6.2),  of  the  case  m  =  \  =  \N  iov  N  =  16, 24, 36, 62. 
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